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We report J/-0 spectra for transverse momenta p T > 5 GeV/c at mid-rapidity in p+p and Au+Au 
collisions at -v/s NN = 200 GeV. The inclusive J/t/j spectrum and the extracted B-hadron feed-down 
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are compared to models incorporating different production mechanisms. We observe significant 
suppression of the J/ip yields for p T >5 GeV/c in 0-30% Au+Au collisions relative to the p+p yield 
scaled by the number of binary nucleon-nucleon collisions in Au+Au collisions. In 30-60% collisions, 
no such suppression is observed. The level of suppression is consistently less than that of high-p T 
7r and low-p T J/ip. 

PACS numbers: 12.38.Mh, 14.40.Gx, 25.75.Dw, 25.75.Nq 



Ultrarelativistic heavy-ion collisions provide a unique 
environment to study strongly interacting matter at high 
temperature and energy density where a transition from 
the hadronic phase of matter to a new partonic phase, the 
Quark-Gluon Plasma (QGP), takes place. Measurements 
of the in-medium dissociation probability of the different 
quarkonium states are expected to provide an estimate of 
the initial temperature of the system [IHi. The J/i/>(lS) 
is the lightest and most abundantly produced quarko- 
nium state accessible in experiment. However, significant 
decay contributions (>40%) from excited cc states, such 
as the Xc and ip(2S), and from B mesons could compli- 
cate the suppression picture suggested by direct dissoci- 
ation models [f| . In addition, other effects absent in p+p 
collisions are likely to have a significant impact on the 
observed J/ip yields in rclativistic heavy-ion collisions at 
CERN-SPS, BNL-RHIC and CERN-LHC [Mil- These 
effects include cold nuclear matter (CNM) effects such 
as initial state parton scattering, nuclear shadowing and 
nuclear absorption, the combined effects of finite J/ip for- 
mation time and the finite space-time extent of the hot, 
dense volume where the dissociation can occur, and re- 
combination of unassociated c and c in the medium fl2| . 
Most of these effects are expected to decrease with in- 
creasing J / ip p T 13J, LLJ] . It is therefore anticipated that 
J/ip measurements at high-p T provide an important tool 
to decouple some of the mechanisms mentioned above 
and provide a cleaner way to extract the contribution 
from color-screening effects 

EMI- 

STAR's previous 

J / ip measurements are consistent with no suppression for 
p T > 5 GeV/c in Cu+Cu collisions at y/s^ = 200 GeV, 
to within the limited precision of the data [13] • In addi- 
tion, the small system size created in Cu+Cu collisions 
may result in high-p r J/ip formation outside the medium. 
Precise measurements in Au+Au collisions are thus cru- 
cial for a systematic study of J/ip production in the hot, 
dense medium. 

The interpretation of medium-induced J/ip modifica- 
tion requires a good understanding of its production 
mechanisms in p+p collisions, which include direct pro- 
duction via gluon fusion, parton fragmentation, and feed- 
down from higher charmonium states and B-hadron de- 
cays The initial hard interactions required to cre- 
ate the charm quark pairs can be well calculated in per- 
turbative QCD (pQCD). However, the subsequent soft 
processes required to form the J/tp hadron and the J hp 
formation time are theoretically not well understood |6j . 
No model at present fully explains the J/tp observations 



in elementary collisions [fjj. The J/ip spectrum in the 
intermediate and high-p r range, together with the angu- 
lar correlations of a high-p T J / ip and associated charged 
hadrons, may provide additional insights in the underly- 
ing production mechanisms. 

In this letter, we report a measurement of J/tp pro- 
duction for 2 < p T < 14 GeV/c in p+p and Au+Au 
collisions at ^/s NN = 200 GeV. The inclusive J/tp cross 
sections at mid-rapidity (|y| < I) in p+p and Au+Au 
collisions are presented. The nuclear modification fac- 
tor, Raa - the ratio of the yield in Au+Au to that in 
p+p collisions scaled by the number of underlying binary 
nucleon-nucleon collisions - is calculated and compared 
to theoretical calculations. 

The Au+Au data used for this analysis were recorded 
in 2010, and the p+p data in 2009. The minimum bias 
(MB) trigger was defined to be a coincidence of the two 
Vertex Position Detectors (VPD) [Isj . Online trigger 
conditions that utilized a MB trigger condition and two 
thresholds for the energy deposited in any single Barrel 
Electromagnetic Calorimeter (BEMC) [l9| tower, with 
a size of A77 x A<p = 0.05 x 0.05, were used to maxi- 
mize the sampled luminosity. To increase the trigger effi- 
ciency, the p+p data with high BEMC threshold were 
recorded without a MB requirement. The p+p data 
with low BEMC threshold were pre-scaled to keep the 
data rate manageable. The integrated luminosities of 
the data samples used for this analysis are 23.1 pb" 1 and 
1.8 pb" 1 with a transverse energy threshold of Et > 6.0 
and 2.6 GeV, respectively, in p+p collisions and 1.4 nb _1 
with E T > 4.3 GeV in Au+Au collisions. In the Au+Au 
data, the collision centrality is determined by the dis- 
tribution of charged-particle multiplicity within \r)\ <0.5 
and Monte Carlo Glauber calculations [20(. 

In this analysis, J/ip —> e + e~ decays were recon- 
structed using the STAR Time Projection Chamber 
(TPC) and the BEMC with full azimuthal cov- 



erage over the pseudorapidity range \rj\ < 1 17|, l23 |. 
Electron identification (cID) for the BEMC triggered 
tracks was achieved by measuring the ionization energy 
loss (dE/dx) and track momentum from the TPC, as 
well as the energy deposition in the BEMC. In addi- 
tion, the shower profile in the barrel shower maximum 
detector (BSMD) [l9[ was used in Au+Au collisions 
to further suppress hadron contamination. At moder- 
ate p T (1 < p T < 3 GeV/c), TPC dE/dx provides 
elD with reasonable efficiency and purity. At low p T 
(0.2 < p T < 3 GeV/c), the elD significantly benefits 
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FIG. 2: (Color online) J/tp p T distributions in Au+Au col- 
lisions with different centralities at ^/s NN = 200 GeV. For 
clarity, the data and curves have been scaled as indicated in 
the legends. The PHENIX results are reported in Q. The 
curves are model fits described in the text. 



FIG. 1: (Color online.) (a) The unlike-sign e + e~ invari- 
ant mass distribution from same-event pairs (filled circles) 
and mixed-event pairs (continuous curve) and their difference 
(open circles) in Au+Au collisions at ^/s NN =200 GeV. (b) 
The invariant J/tp cross section versus p T in p+p collisions at 
y/s = 200 GeV. The vertical bars and boxes depict the statis- 
tical and systematic uncertainties, respectively. Also shown high S/B ratio. About 400 J/i/j with p T > 4 GeV/c and 



are results published by STAR [13] and PHENIX [21j]. The 
curves are theoretical calculations described in the text, (c) 
The fraction of B — > J/tp over the inclusive J/tp yield in 
p+p collisions. The FONLL+CEM model calculation is also 
shown. 



from a recently installed large area time-of- flight (TOF) 



detector covering \t]\ < 0.9 [24426j . The complete TOF 



detector was available for the 2010 Au+Au run, whereas 
72% was available for the p+p data collected in 2009. 

The J/tp signal was extracted by subtracting from 
the unlikc-sign ee invariant mass spectrum the random 
combinatorial background that was reproduced by the 
like-sign spectrum in p+p collisions and unlike-sign spec- 
trum from mixed-events in Au+Au collisions [27] . Fig- 
ure [IJa) shows the invariant mass distribution before 
and after the combinatorial background subtraction in 
Au+Au collisions at y/s^ = 200 GeV. The J/ip raw 
yields were obtained from a mass window of 2.7 < Mf^ v < 
3.3 — 3.4 GeV/c 2 in p+p collisions depending on the J/tp 
p T , and 2.9 < M?£ v < 3.2 GeV/c 2 in Au+Au collisions. 
The yields were corrected for « 10% radiation losses that 
cause some of the decay daughters to be reconstructed 
with Mf^ v outside the above mass ranges. The total J/tp 
yield was sa 1100 (p T > 2 GeV/c) in p+p collisions. It 
was w 1000 and 300 {p T > 3 GeV/c) in 0-20% and 40-60% 
Au+Au collisions, respectively. The signal to background 
ratio (S/B) was « 4 in p+p collisions and « 1/7 (1/2) 
in 0-20% (40-60%) Au+Au collisions. In the J/tp-hadron 
correlation analysis, the cuts were selected to provide a 



3.0<M?„ e u <3.2 GeV/c 2 were observed with a S/B ratio 
of 22/1 in p+p collisions [Hj|. 

Acceptance and efficiency corrections were studied us- 
ing Monte Carlo GEANT simulations 0. The system- 
atic uncertainties are due to kinematic cuts (< 12%), sig- 
nal extraction (including the contribution from correlated 
background, 2-26%), momentum resolution (<3%), and 
efficiency (7.5%). The normalization uncertainty for the 
cross section in p+p collisions is 8.1% (2j|. The corre- 
lation of the systematic uncertainties in the spectra was 
taken into account in evaluating the uncertainty on Raa- 
FigurcQJb) shows the J/tp invariant cross section times 
the braching ratio (Bu) |30| as a function of p T for p+p 
collisions at y/s = 200 GeV. The new results are consis- 
tent with those previously published by STAR [l?} and 
PHENIX [HI. The rapidity window for PHENIX results 
is \y\< 0.35. The dashed curves depict ncxt-to-leading or- 
der (NLO) theoretical Non-Relativistic QCD (NRQCD) 
calculations from color-octet (CO) and color-singlet (CS) 
transitions [31( for prompt J / tp production in p+p colli- 
sions. The CS+CO calculations match the p T spectra for 
p T > 4 GeV/c to within the uncertainties. The continu- 
ous curve shows the calculation from the color evapora- 
tion model (CEM) for prompt J/tp (3^]. It describes the 
p T spectra reasonably well. The bottom panel shows the 
ratios of the data and theory calculations to the CEM 
calculation. A model based on NNLO* CS 13311 predicts 
a too steep p T dependence, as discussed in [17J. 

The relative contribution of £?-hadron feed-down to 
the inclusive J/tp yield was obtained in the same way as 
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FIG. 3: (Color online.) J/tp Raa versus p T for several cen- 
trality bins for Au+Au collisions at ^/s NN = 200 GeV. The 
statistical (systematic) uncertainties are shown with vertical 
bars (open boxes). The filled boxes about unity on the right 
show the size of the normalization uncertainty. PHENIX low- 
p T J/tp results [D] and STAR high-p T results in Cu+Cu col- 
lisions [l7j are shown for comparison. The curves are the 
predictions by Model I (Liu et al.) [ll| and Model II (Zhao et 
al.) 0. 



FIG. 4: (Color online) Raa versus N par t for high-p r J/tp, 
low -£x, J H from PHENIX @], and high-p T tt* from STAR 
EEL Eg ■ The statistical (systematic) uncertainties are shown 
in vertical bars (boxes). The shaded green band about unity 
shows the systematic uncertainties from Nun and the box 
about unity on the right shows the -Raa normalization uncer- 
tainty from the statistical and global systematic uncertainties 
of the p+p reference data. 



m [17J, |23|. We note that this method is data-driven, 
although it relies also on the validity of PYTHIA's 
3J, |35j modeling of the near-side associated hadron 
distributions. The effect is found to be 10-25% for 
4 < p T < 12 GeV/c as shown in Fig. (He). Within err- 
ors our data are consistent with the Fixed Order plus 
Next-to-Leading Logarithms (FONLL) plus CEM pre- 
diction 3^, 37 [ indicated by the curve and uncertainty 
band. More precise measurements using displaced vertex 
techniques 38, 3i| similar to those employed by CDF in 
p + p collisions at y/s = 1.96 TeV and by ATLAS and 
CMS in p+p collisions at y/s = 7 TeV [40^42] are needed 



37 



to quantify the anticipated energy dependence [36 

The measured J/tp p T spectra in Au+Au collisions for 
different centralities are shown in Fig. and are consis- 
tent with the low-p T J/tjj results from PHENIX in the 
region of overlap in p T . The continuous curves depict fits 
based on the Tsallis statistics Blast-wave (TBW) model 
to the combined STAR and PHENIX data with the ra- 
dial flow velocity /3 fixed to zero (13] • The fits reproduce 
the data reasonably well. Under the assumption that the 
J/tp flows like light hadrons [H, Hj] , the TBW calcula- 
tions shown as dashed curves underpredict the yields at 
low p T . This could be due to a small (or zero) radial 
flow or a significant contribution from charm quark re- 
combination that would enhance the yield at low p T , or 
both. 



Figure [3] shows the J/tp Raa versus p T for different 
centrality bins. The STAR Cu+Cu J/tp results [l7| are 
also shown. For p T > 5 GeV/c, J/tp Raa m the 40- 
60% centrality bin is consistent with unity and with our 
previous measurement in Cu+Cu collisions with a simi- 
lar average number of participants (N var t)- The curves 
show two theoretical calculations 13, [lj] describing the 
data reasonably well. These calculations include contri- 
butions from prompt production and statistical charm 
quark regeneration. The suppression of the prompt J/tp 
component in the model calculations is mainly due to 
the color-screening effect. The model from Zhao et al. 
(Model II) [14j also includes the J/tp formation-time ef- 
fect and the i?-hadron feed-down contribution. 

For p T > 5 GeV/c, J/tp production follows the scal- 
ing of the cross section, = g{xr) / {\fs) n , with 
xt = 2vt/\/s [47l [ill ] observed in p+p collisions at t/s 
— 200 GeV [17J. This indicates that soft processes do 
not alter high-p r J/tp production. In this p T region, the 
CNM and cc recombination effects are expected to be 
negligible in heavy- ion collisions [H, 14|, and AdS/CFT 
predicts suppression in the QGP for direct J/tp [l5j. 

We present high-p T (p T >5 GeV/c) Raa as a function 
of Npart in Fig. |4] No significant suppression of high-p r 
J/tp production is observed in mid-central to peripheral 
collisions (30-60%, N part < 140). In central collisions 
(0-30%, N par t ^ 140), high-p r J/tp are significantly sup- 
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pressed. The Raa of low-p r (0 < p T < 5 GeV/c) J/ip 
measured by PHENIX Q and higW (p T > 5 GeV/c) 
charged pions measured by STAR [45|, |46[ are shown for 
comparison [1, HBl [Zfjj]. The high-p T J/i/> Raa is sys- 
tematically higher. The predictions of high-p T J /if) Raa 
from Model I (Liu et al.) |H| and Model II (Zhao et 
al.) [IH are shown as continuous and dashed curves, re- 
spectively. Model I describes our data reasonably well. 
Model II underpredicts J/ip Raa at N part > 70. 

The high-p T J/ip Raa versus centrality is different 
from that of high-p T pions. This is expected from dif- 
ferences in their production. Dissociation is considered 
to be the dominant mechanism in the case of J/ip pro- 
duction, and induced gluon radiation in the case of pion 
production. For p T > 5 GeV/c, the recombination and 
initial parton scattering effects are expected to be neg- 
ligible [13, LL4l • The observed J/ip Raa dependence on 
p T and system size might be due to the interplay of for- 
mation time, color screening, sequential suppression, and 
parton distribution functions in heavy nuclei [l2j . Future 
J / ip measurements in d+ Au collisions over a broad range 
in p T are needed to understand the CNM effects at RHIC. 
The i?-hadron feed-down contribution is assumed to be 
the same in Au+Au and p+p collisions. Based on our 
measurement, Raa for prompt J/ip with p T > 5 GeV/c 
in the most central collisions will be 0.80 if B-hadrons 
are completely suppressed and 0.55 if B-hadrons are un- 
supprcssed. 

In summary, we report measurements of J/ip produc- 
tion in y/s^f = 200 GeV p+p and Au+Au collisions 
for p T > 2 — 3 GeV/c at RHIC. The p T spectrum in 
p+p collisions is compared to various theoretical cal- 
culations. Currently, only the CEM model and NLO 
CS+CO calculation describe our data. Based on the 
measurement of azimuthal correlations between high-p r 
J/ip and charged hadrons we estimate the fraction of J / ip 
from 5-hadron decay to be 10-25% in the p T range of 
4-12 GeV/c in p+p collisions. The nuclear modification 
factor Raa in Au+Au increases from low to high p T . For 
p T >5 GeV/c, J/ip Raa is consistent with no suppression 
from mid-central to peripheral collisions (30-60% central- 
ity), and significantly smaller than unity in the most cen- 
tral Au+Au collisions. The results on Raa versus p T and 
Npart provide new insight in the study of color screening 
features for charmonium. 
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